Exact results on the two-particle Green's function of a Bose-Einstein condensate 
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Starting from the Dyson-Beliaev and generalized Gross-Pitaevskii equations with an extra nonlo- 
cal potential, we derive an exact expression of the two-particle Green's function 1C for an interacting 
Bose-Einstein condensate in terms of unambiguously defined self-energies and vertices. The formula 
can be a convenient basis for approximate calculations of JC. It also tells us that poles of IC are 
not shared with (i.e. shifted from) those of the single-particle Green's function, contrary to the 
conclusion of previous studies. 
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The realization 1 of Bose-Einstein condensation (BEC) 
with an atomic gas in 1995 has revived intense the- 
oretical interests on interacting condensed Bose sys- 
tems. One of their unique features is that the gapless 
Nambu-Goldstone boson^ of the broken U(l) symme- 
try, i.e. the Bogoliubov mode^ emerges as a pole of the 
single-particle Green's function G to dominate thermo- 
dynamic properties. It also seems to have been widely 
accepted that poles of G are shared with those of the two- 
particle Green's function jC, as first claimed by Gavoret 
and Nozieres^ in 1964 and reproduced by the dielectric 
formalism^— These theories have provided a support to 
utilize G for describing collective modes of condensed 
atomic gases. Indeed, the sharing of common poles be- 
tween G and jC has been regarded as one of the most 
spectacular features of condensed Bose systems. 

However, the theory by Gavoret and Nozieres^ is based 
on an analysis of the structures of simple perturbation ex- 
pansions performed separately for G and JC. Thus, it may 
suffer from ambiguity as to how to define self-energies and 
vertices in the presence of an "improper" interaction hav- 
ing only a single quasiparticle channel inherent in BEC. 
Since BEC is a prototype of broken symmetry, it will be 
well worth reinvestigating the fundamental issue with a 
different method and viewpoint. 

As is well known in normal systems^— a two-particle 
Green's function can be generated from a single-particle 
Green's function by a functional differentiation with re- 
spect to an additional potential. This method enables us 
to derive a formally exact expression of the two-particle 
Green's function in terms of unambiguously defined self- 
energies and vertices. Moreover, it can be used in practi- 
cal calculations of the two-particle Green's function with 
Baym's ^-derivable approximation.-lfl The approximation 
has a great advantage that the whole series of thermo- 
dynamic, single-particle, and two-particle properties can 
be discussed in a unified way based on a single functional 

even beyond equilibrium. 11 

We here apply the functional-differentiation method 
to an interacting Bose-Einstein condensate to obtain an 
exact expression of JC. The formula can also be used 
for practical calculations of JC with the self-consistent 
^-derivative approximation of condensed Bose systems 



developed recently.— Our derivation is based solely on 
rigorous results of the previous paperJ^ It will thereby 
be shown that poles of JC are not shared with those of G, 
contrary to the previous conclusion^— Unlike the previ- 
ous studies for homogeneous systems using the momen- 
tum conservation^ — our formulation will be carried out 
in the coordinate space so that it is applicable to trapped 
atomic gases. 

We consider identical Bose particles with mass m and 
spin described by the Hamiltonian: 



H = J dW(ri)^i<A(ri) + ~ J £ 

XiI) ] (y 2 )V{y 1 - V 2 )tp{Y2)lp{v 1 ). 



n / d 3 r 2 ^ t (ri) 
(1) 

Here tp^ and ip are field operators satisfying the Bose 
commutation relations, K\ = —h 2 'Vi/2m — [i with /i the 
chemical potential, and V is the interaction potential. 
Though dropped here, the effect of a trap potential can 
be included easily in K \ . Let us introduce the Heisenberg 
representations of the field operators by 



^i(l) 



V>(ri)e 



-TiH 



^(1) = e^VMe -71 * 



(2) 

with 1 ee (ri,ri), where < n < T" 1 with T the tem- 
perature in units of h = = 1. The operators V'i(l) 
and ^(l) were denoted previously^ by tp(l) and V'(l), 
respectively. We next express tpi(l) as a sum of the con- 
densate wave function \&i(l) = (tpi(l)) and the quasipar- 
ticle field fa(l) as 



= + &(l) (i = l,2), 



(3) 



with (• • • ) the grand-canonical average in terms of H. 
Note: (i) (fa(l)) = by definition; and (ii) ^i(l) = 
^2(1) = ^( r i) m equilibrium with the superscript * sig- 
nifying complex conjugate. Using fa, we introduce our 
Matsubara Green's function in the 2x2 Nambu space 



G«(l,2) = -(T T i (l)</)3- i (2)>(-l^- 1 , 



(4) 



where T T denotes the "time" -ordering operator^ They 
satisfy^ 



0^(1,2) = (-l) i +'- 1 G 8 _ J ,3-i(2 ) l) 
= (-l) i +^G^(r 2 r 1 ,r 1 r 2 ). 



(5) 
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Let us recapitulate exact results on the matrix G — 
(Gjj) and the vector * = [\f r i\E r 2 ] T ; see Sec. II of Ref. 
\12\ for details. First of all, they obey the Dyson-Beliaev 
equation and the generalized Gross-Pitaevskii equation 
(or generalized Hugenholtz-Pines relation) given by 



G^ 1 (l,3)G(3,2) = a <5(l,2), 



G- 1 (l,2)a 3 V(2) = 0, 



(6a) 



(6b) 



respectively. Here summations over barred arguments 
are implied, &o an d 63 denote the 2x2 unit matrix 
and the third Pauli matrix, respectively, 6(1, 2) = 8(t\ — 
t 2 )<5(i"i — r 2 ), and G -1 is defined by 

G~ l (l,2) = (-ao-^-asKx) 5(1,2) -£(1,2), (7) 

with E the self-energy matrix. We point out that the 
first component of Eq. f|6b[) in equilibrium is written ex- 
plicitlyas-Mtr!) =Eii(l,2)*(r 2 )-E 12 (l,2)**(r 2 ). 
By approximating En(l,2) = 2g8(l, 2)|*(ri)| 2 and 
E 12 (l, 2) = gS(l, 2)[*( ri )] 2 for V(v x - r 2 ) = g5(v 1 - r 2 ), 
it reduces to the standard Gross-Pitaevskii equation.— ~— 
Setting \I/ — » \/n and K — > —ji with no the con- 
densate density in the same equation, we also ob- 
tain the Hugenholtz-Pines relation for the homogeneous 
system i 12 i 17 

It has been showr*i£ that the elements of E satisfy the 
same relations as Eq. JS]). Moreover, all of them can be 
obtained from a single functional $ = $[G, F, F, ^ 2 ] 
as Eq. (21a) of Ref. 01 with G = G n , F = G 12 , and 
F = — G 2 i. Using Eq. (0, we here write every G in $ 
as G(l, 2) = [G u (l, 2) - G 22 (2, l)]/2. Then the relevant 
relations can be put into the single expression: 



Sy(l,2) = - 



2 <S$ 



(8a) 



T5G 3l {2,iy 

The functional $ also satisfies Eq. (21b) of Ref. [l2l i.e 
1 5<P 



Ts*7W) =J * ilM - ir1 " 3 ®- (8b) 

With these preliminaries, we now study the two- 
particle Green's function: 

£y,iw(12,34) 

= (T T tpi(l)ipk (3)^3_, (4)^3-, (2)) (-iy +L 

-(T r ^(l)^-,(2))(T T ^(3)^3-K4))(-l) J+ '. (9) 

Collective modes correspond to the poles of this Green's 
function. To derive the equation for K,, we follow a 
standard procedure to produce the two-particle Green's 
function from G££ Let us add an extra perturbation de- 
scribed by the S matrix: 



S((3) = T T exp 



•^?(l)V3-j(2)(- 



U Ti (2, 1) 



with /3 = T -1 . The full Matsubara Green's function 
in the presence of the nonlocal potential U = (£/y) is 
defined byi^ 

(T rt S(/3)^(l)^ 3 - 3 (2)) 

^•^ 2 ) = jsm (-^ 

= -(T T 5(/3)^(1)^3-,(2)) C (-1) J - 1 , (Ha) 

where the subscript c denotes contribution of those Feyn- 
man diagrams connected with ipi(l) and/or ^3__y(2). 
Noting that there may be the finite average \Et(l) = 
(T T S((3)i/ji(l)) c , we can transform Eq. (|llal) into 

^•(l,2) = G <i (l J 2)-* < (l)* 3 _ i (2)(-iy- 1 , (lib) 

with G y -(1,2) = -(T r 5(/3)0 i (l)03-,(2)) c (-l)^ 1 ; this 
quantity reduces to Eq. Q as U — > 0. It may be seen 
easily that two-particle Green's function (|9]) is obtained 
from Eq. (jllaj) by 



/Cy, w (12,34) = 2 



^■(1,2) 
5U lk (4,3Y 



(12a) 



where the limit U — > is implied after the differentiation; 
we will use this convention below. A substitution of Eq. 
(ITTbf into Eq. (fl"2"a|) yields 



£^,(12,34) = 2 



+ 



5G l3 (l,2) 
6U lk (A,3) 

6U lk (4,3) 



- 2 



6U lk (A,3) 



*a-i(2) 



(12b) 



Equation (|12b|) tells us that we only need to know the 
linear responses of G and ^ to U for writing JC down 
explicitly. 

To carry it out, we start from Eq. (|6|). Differentia- 
tions of Gy(l,2) = -(T r S(/3)^(l) ( fe- J (2))c(-l)''- 1 and 
ee (T T S(/3)ipi(l)) c with respect to n tell us&£ that 
perturbation, Eq. (flU)) . adds to the right-hand side of Eq. 
© an extra term -U'(l,2) with 



v ,, 12) = ^■(L2) + (-l)^- 1 £/ 3 ^3-z(2,l) ^ 



(13) 



Varying [/—>[/ + <5J7 and subsequently setting [7 = in 
resultant Eq. (|6]), we obtain the first-order equations, 

G- 1 (1,3)5G(3,2)= [6U' (1,3) + 6±(1,3)]G(3,2), (14a) 

G- 1 (l,2)a 3 ^(2)= [6U'(l,2) + 6£(l,2)]a 3 $(2). (14b) 

At this stage, it is convenient to introduce the following 
quantities: 



(4) , , 1 £S;,-(1,2) 1 

T% (12,34) ee 



6 2 $ 



(10) 



2^(4,3) T SGji(2, l)5G/fe(4, 3) 

(15a) 
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)8C )©- -@C m Using r ^ (12, 34) = ( - 1)i+j_lr 3 4 Vi,«( 21 ' 34 ) 



FIG. 1: Irreducible vertices of Eq. (fT5j) . 



ra(i2,3).i(-i)^ 5E M(M) 



<5* fe (3) 

(4) / 
ij,3-k,V 



2(-l)*+V 4 > (12, 34)^(4), (15b) 



fg fc (l,23) = 2(-l)'- 1 r ( g fe (14,23)* r (4) 

= (-i) 4 rf fc ) , 3 - 4 (23,i) I 



(15c) 



rg ) (l,2)^2(-l) S - 1 r^.(13,2)vI/ s (3), (15d) 

where Eq. ([5]) has been used to derive the second ex- 
pression of r^ 4,3 '. These are "irreducible" vertices of 
our condensed Bose system, as seen below, and can be 
expressed diagrammatically as Fig. [1] It follows from 
Eq. <j5j) and $* = $ that they satisfy various sym- 



metry relations, e.g., r,^ fc; (12,34) 



034,12) 



(-l)* + *'- x r£Zj,3_ i>w (21,34). The quantities and 
p(3) correspond to / and J of Gavoret and Nozieres^ 
respectively. Our definitions may be advantageous over 
theirs because the vertices can be obtained explicitly 
from a single functional $ with clear relations among 
them. 

Using r' 4,3 ) above, we can express ST, in Eq. (fT4"|) as 

SZijQ.,2) = ^^(12,43)^(3,4) 

+2r£(12,3)(-l) S - 1 <^(3> (16) 

It enables us to transform Eq. ([T4"]) into a closed set of 
equations for SG and <5\E f . Indeed, multiplying Eq. (fT4")) 
by G from the left, substituting Eq. (fT6) , and using Eqs. 
(IT5c|) and (fl"5dl) . we obtain 

«?y(l,2) 

= G ir (l,4)G £i (3,2)^(4,3) 

-2G ir (l, 4)G Si (3, 2)rgJ„ (43, 65)^^(5, 6) 

+2G 4r (l,4)G £j (3,2)rg ?fi (43,5)(-l)' ft - 1 ( 5*„- l (5), 

(17a) 



ij,fe2 \ ' / V / 3— Jj3— i,kl * 

we can also transform G^(l, 4)G^ ; (3, 2)r|| ) mn (43, 56) 

(-l)^- 1 G ii3 _ £ (l I 3)G 3 _ rii (4,2)rg mn (43,56). 

To provide Eq. (fl~7)) with a compact expression, let us 
introduce the vectors 8G and <5C7 by 

(12^1^ = ^(1,2), (12^1^ = ^(1,2), (18) 

together with the matrices K, T (4) , x ( °\ L E (3) , f (3) , 
^ 3 M (3) , and f( 2 ) by 

(12«|£|43 Jfc ) =/Cy,w(12,34), (19a) 



(12 y |r( 4 )|43 ;fc )^r^(12,34), (19b) 



(12 4 ,|x ( °)|43 ifc ) ee G«(l,4)G fcj -(3,2) + (-l)^'" 1 

xG i) 3- fc (l,3)G 3 -i J -(4,2), (19c) 

(12y|l|43 a ) =6u5 kj S{l, 4)6(3,2), (19d) 



<i2^|r( 3 )i3 ft )=irgici2,3), 



(l i |f^|32 fcj )=r^(l,23) ) 



(19e) 



(19f) 



(12 y |*( 3 )|3 fe ) ee (-1)^^(1)4,3-^(3,2) 

+5 w *(3,l)*3-i(2)] J (19g) 



(l i |i (3) |32 fci ) ee (-I)'" 1 [^(1,3)^(2) 

+* j)3 _ J -*(l s 2)*3-fc(3)], (19h) 

(1^12^ = ^(1,2). (19i) 

The quantity x'" '' describes independent propagation of 
two particles. 

Using Eqs. (fT8)) and (fT9|) and noting the comments be- 
low Eq. (|17bj) . we can express Eq. (17]) as JG = |x (0) <5f7- 

X (0) r (4) 5G + x (0) r (3) <7 3 (5* and <7 3 6# = ±G*i (3 W - 

GF (3) 5G+Gf ( 2 )(T3(5*. They are further transformed into 



(-lr^i) = ^(Mx-iy-^^M 

-%(l,2)f^(2, 43)<JG Er (3, 4) 

+Gy(l,2)r|J(2,3)(-l)*- 1 «* s (3). 

(17b) 

Note G jr (l,4)G^-(3,2)^(4,3) = |[G a -(l,4)G Si (3,2) + 
(-l) fe +'- 1 G 4! 3_ s (l,3)G 3 _ rj (4,2)]W s (4,3) from Eq. 



£G = ^x (4) <^ + x (4) r (3) <r 3 (5*, (20a) 



(20b) 



where x & n d X^ "* are defined by 

x (4 ) ee (i + x (0) r (4) rV 0) = ix^- 1 + r^y\ (2ia) 
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FIG. 2: Diagrammatic representation of Eq. (|21[) . Every long 
straight line in the second equation denotes G. 



x (2) = (i- Gfwy'c = {g- 1 - f < 2 >)~\ 

It is also convenient to introduce 



(21b) 



(21c) 



^-fW+f^xW 35 ) -1 , (21d) 



where the superscripts q and c denote "quasiparticle" 
and "condensate," respectively. Figure 2 expresses Eqs. 
(|21a[) - (|21d|) diagrammatically. Now, we can write down 
the solution to Eq. (|2T)]) in terms of x^ and x as 

SG = ^x (q) (I + r^x^± i3) )SU, (22a) 

^ = ^ 3 x (c) (i (3) -£ (3 V 4) )^- (22b) 



Let us substitute Eq. (f22|) into Eq. (|12b[) and make use 
of Eq. (HU) as well as x (q) r (3) x (2) = x (4) r (3) x (c) in Eq. 
0}- We thereby obtain JC defined by Eq. (|19aj) as 



ic = x (q) + x (4) r (3) x (c) * (3) + * (3) x (c) f (3) x (4) 



X + X 

^(3)^(c)J,^ 



(23) 



This expression clearly tells us that collective modes are 
determined as poles of x an d X^- Note in this con- 
text that poles of x^ m Eq- (|23|) are cancelled by those 
of x'' 4 '' in the denominator of x^ c \ as seen from Eq. 
(f2"Td|) . It also follows from Eq. ([2l"dl) that the poles of x (c) 
are generally not identical to those of the single-particle 
Green's function G due to the additional contribution 
f^ 2 ' — r^x^r 1 - 3 '', in contradiction to the statement by 
Gavoret and Nozieres^ This point will be discussed in 
more detail below. 

Equation (gH]) with Eqs. (dHJ), (0, and (ED is tne main 
result of the present paper. The expression is formally 



exact, clarifies the structure of the two-particle Green's 
function IC in terms of unambiguously defined vertices, 
and enables us to carry out practical calculations of IC for 
a given approximate $ on the same footing as thermo- 
dynamic and single-particle properties^ The last point 
may be regarded as a definite advantage of the present 
formalism over the dielectric one£ — 

Equation (|23[) in the coordinate representation can be 
used to investigate two-particle correlations of general in- 
homogeneous systems, including homogeneous ones. For 
the latter cases, however, it is far more convenient to 
adopt the "energy" -momentum representation. To be 
specific, vertices (fT5|) in those cases can be expanded as 



pp'q 

s ip'-f 3 -i(p'+q)-f4 



rg>2,34) = Y.^&p'^^^^ 

(24a) 



xe ' 



a i(p+q)-ri-ip-r 2 -iq-f 3 



, (24b) 



p'q 



iq-ri+ip-r2—i{p+q)-r 3 



rg ) (i,2) = ^rg ) (^e^ 



(ri-f 2 ) 



■ r3 , (24c) 



(24d) 



where fi = (ri,— ti), p = (p,e n ) with e n = 2mtT 
(n = 0,±1,---), and the summation over p denotes 
TY, n fd 3 p/(2ir) 3 . Other quantities in Eq. $f§§ can be 
expanded similarly. The Fourier coefficients of Eqs. (|19c[) . 
(|19dl) . ([TOgf, and (I19h|) are thereby obtained as 



xfiuiP.p'.q) = 5p l pG ll {p + q)G k ,(p) + (-1) 

x<V -p*-q-Gi,3-k(p + <z)G 3 _ ij (p),(25a) 

hj,ki(P^P'^ = &uhj5p>p, (25b) 
*;&(££) = (-V* +k >/nS(h,3-j6 1 r,-t+6H8 f tf), (25c) 
*i%(p, <t) = (-i^^Mp-g + <v, A- ■;) > (25d) 

respectively, where 5p~'p~ = (2ir) 3 T~ 1 5(p' — p)S n ' n and no 
denotes the condensate density. It then follows that Eqs. 
(f2"Tj) and (|2"3"j) hold as they are in terms of the Fourier 
coefficients. For example, Eq. (|21c[) can be written ex- 



k+l-l 



plicitly as an integral equation for Xw fcj(pj p' '■> Q) as 



X {q) (P,P',Q) = X (4) (P,P',<1) ~J2x {4) (p,Pi,q)L {3) (pi,q) 

P1P2 

X X (2) (q)t 3) (P2, q)X (q) (P2 , p', q) , (26) 
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where X^ q \ etc., are now matrices only in terms of the 
Nambu indices i, j, ■ ■ ■ , which may be defined explicitly 
as Eq. (|19j) without space-time arguments. 

We now compare Eqs. (f2"Tj) and ([25)1 with the results 
for the two-particle Green's function JC by Gavoret and 
Nozieres. 4 Apparently, they found the same structure for 
JC as Eq. ([25)1 above. They subsequently identified the 

quantity corresponding to S+f ( 2 ) — X^T^ in \^ of 
Eq. (I21dp with the single-particle self-energy as Eq. (3.4) 
of their paper, where M and JGGP on the right-hand 

side correspond to £ + f^ 2 ' and — £ ^X^E^j respec- 
tively. However, they did not provide detailed reason- 
ing to the crucial statement. In this context, we would 
like to point out that their analysis of JC was carried 
out separately from that of G by only investigating its 
diagrammatic structure in the simple perturbation ex- 
pansion, where r^ 2 -*, for example, may be mistaken eas- 
ily for a part of the single-particle self-energy, as seen 
from the second diagram of Fig. 2. This may be the rea- 
son why they concluded erroneously that the quantity 

corresponding to t + f (2) - r (3) x (4) E (3) is the single- 
particle self-energy. In contrast, our investigation of JC 
has been performed on the basis of Eq. © for G and 'S, 
where the self-energy 2 is defined unambiguously at the 
single-particle level. It is thereby shown that the term 



f ( 2 ) _ f ^(4)p(3) ghc-uid be regarded as additional con- 
tribution distinct from the single-particle self-energy. 

Thus, of fundamental importance will be to clarify how 
the extra "self-energy" f< 2 ) - f (3) X (4) E (3) in X (c) shifts 
its poles from those of G. In the weak-coupling limit, 
we can show = and ^(1,2) = -r^*(l,2) = 
-2VTri - r 2 )*(ri)*(r 2 ) by using Eqs. (25) and (26) of 
Ref.Qjand Eq. JTjjJ) above. Combined with £i 2 (l,2) = 
V(ri - r 2 )[#(ri)tf(r 2 ) - (0(1)0(2))] from the lowest- 
order gapless <l>-derivable approximation^ we thereby 
obtain Si 2 (l, 2) + r 12 (l, 2) = V(r 1 - r 2 )[-*(n)*(r 2 ) - 
(0(1)0(2))]. Thus, at the mean-field level, f( 2 ) merely 
changes the sign of the condensate (i.e., dominant) con- 
tribution to the off-diagonal self-energy. It hence fol- 
lows that, to the leading-order in the interaction, poles 
of (G^ 1 — r' 2 ))" 1 are the same as those of G. How- 
ever, they are not exactly identical due to the presence 
of (0(1)0(2)). Beyond the weak-coupling regime where 

the polarization contribution — f X^E a ^ so becomes 
relevant in x^ c \ therefore, it is reasonable to expect that 
poles of x an d G are generally different. Further in- 
vestigation needs to be carried out on the similarity or 
difference between the single-particle and collective exci- 
tations. 
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